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A seemingly natural mechanism is proposed, that could stop the gravitational col- 
lapse of a very massive body. Without needing to change the concept of the collapsing 



process itself, that is, without invoking thin layers nor resorting to asymptoticity (as 
has been usually done in the literature), it is proven that a model can be built in which 



$h , the quantum vacuum is able to produce a negative stress that may stop the collapse of 

the black hole, reaching a final state of the spacetime structure that is a static de Sitter 



model. The solution is found by looking into a generic family of spacetimes: that of 
maximal spherically symmetric ones expanded by a geodesic radial null one-form from 
flat spacetime. They are called here GRNSS spaces, and are proven to constitute a 
distinguished family of Kerr-Schild metrics. The models considered previously in the 
literature are easily recovered in this approach, which yields, moreover, an infinite set 
of possible candidates for the interior of the black hole. First steps towards their semi- 
classical quantization are undertaken. It is shown that the quantization protocol may 
be here more easily carried out than within conformal field theory. 
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1 Introduction 



Soon after the appearance of General relativity, K. Schwarzschild found a solution of Ein- 
stein's equations [jy that was useful to describe both the exterior and the interior spacetime 
metric created by a spherically symmetric body in the case of uniform energy density. This 
very early model already showed a new feature of General Relativity that still concerns the- 
oretical physicists. According to General Relativity there exists a limit involving the mass 
and the size of a star above which, in accordance with classical matter properties, it can- 
not exist and the appearance of a singularity is unavoidable @-||. In the late 60's, S.W. 
Hawking and R. Penrose gave a general framework to analyze gravitational collapse, and 
established their now well-known Singularity Theorems (see e.g. 0). This seemed to settle 
down the issue, but it has been suggested since then, for some decades, that other effects, 
most likely those coming from quantum mechanics ( see e.g. H7|~IP-0|) might alter those pre- 



dictions in a significant way (see e.g. ||11|| - -||17||). Nowadays, the issue of the non-avoidance 
of singularities of space-time is still far from being resolved. However, there are more and 
more indications that, indeed, quantum effects may become most principal agents and that 
a regularized object instead of a true singularity could be formed after the collapse of a very 
massive object. 

As everybody knows, Einstein's equations, or those corresponding to the low energy limits 
of the usual candidates for suitable quantum theories of gravitation, have two distinguished 
parts. One corresponds to the geometry and the other to the stress-energy content. It is 
clear that, for any solution to exist, the perfect matching of both pieces is essential. Due 
to the difficulty in getting (regularising) the contribution from the stress-energy tensor, it is 
usually assumed that the spacetime can still be described classically without much loss. In 
this sense, our work (divided in two papers) gives a general solution for the geometry that 
may be expected to be born after the collapse of a non-rotating object takes place, following 
both old and recent ideas on the type of core that is formed and on the precise way this 
happens (see e.g. fCB|-PB|). The case of rotating objects will be considered elsewhere [ |53|| . 

The case dealt with in the present paper is only the first step of the whole approach, 
corresponding to the resolution of the left hand side of the aforementioned equations. Its 
purpose is to begin the process of building up a way to find, within the distinguished (and 
quite generic) family of solutions that will be obtained in the first part of the paper, some 
compatible quantum fields that, once regularized, may yield the same result as, at least, 
a particular energy-momentum tensor inside this family. As we shall see, the variety of 
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possibilities that comes out, in addition to those of some previous attempts by different 
authors (also recovered here), and the good compliance with basic quantum requirements, 
compels us to believe that this goal can be attained and, in particular too, that black hole 
singularities might in the end be indeed removed by quantum effects. 

The present paper deals with the most plausible models, whereas in a second one (referred 
to as paper II in what follows) it will be proven that they constitute indeed the general 
solution for the stress-energy that was assumed to exist. It is seen also, that isotropization 
takes place far from the regularization scale. We will also deal there with the possible 
extension of the formalism to other candidates, e.g. stringy black holes. 

This paper is two-fold. On one side, as mentioned elsewhere, we show in detail how the 
obtaining of regularized sources for the interiors of non-rotating black holes can be carried 
out, in fact, without changing the initial concept of the collapsing process itself, i.e. without 
having to invoke thin layers nor very detailed asymptotic behaviors to get the expected 
results only there. For the general case of an arbitrary equation of state, we refer the reader 
to [p9| , where some explicit solutions are constructed, yet with a different type of energy- 
momentum tensor than the one we consider here. In Sect. we start with a generic family of 
spacetimes, which is wide enough to include all those that will be relevant for our purposes. 
In Sect. |3|, we write down the matching conditions in order to couple any pair of spacetimes 
of the family under study through an arbitrary hypersurface. We get the result that the 
physical solution can only be a three-sphere of constant "radius". In Sect. f|, we briefly 
review the impossibility of matching Schwarzschild and de Sitter spacetimes directly. In 
Sect. |5|, we specify the matching conditions for a "classical" model of a Schwarzschild- de 
Sitter exterior. And in Sect. || we write the correct solution for the matching of a de Sitter 
interior, as the core of the object, with an exterior which fits with current observational data. 
One easily checks that the de Sitter core is the least convergent regular solution, so that it 
is the limiting solution between singularity- free and singular cores. In Sect. [7], we prove that 
the transition expected to take place is accomplished, provided we demand regularity of the 
collapsed body only. We recover the two known models from our analysis, and show that 



there is indeed an infinite set of possible interiors (see also pi]). In Sect. g[ we carry out 
a numerical analysis for all the candidates of the preceding section. The results show that 
the object is mostly far away from the regularization scale where the notions of space and 
time lose their sense, and therefore it still may be analyzed with some approximate methods 
in order to find out possible quantum sources. Before doing this, in Sect. [| some aspects 
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regarding the horizons and the topological structure of the solutions are briefly considered. 
In Sect. we study the energy conditions in these sources, and show that they can be 
easily related to the properties of the energy density only. In Sect. [II], we start the study 
of the sources that may yield some of the solutions of the general family. In this respect, 
we make the first steps towards their (semi-classical) quantization. Preliminary results show 
that they are indeed good candidates, deserving further analysis, being on the same footing 
as when one deals with a conformal field theory. At the same time, we show that they 
fulfill the conditions for an anomaly origin of their sources. We also indicate how quantum 
field corrections can be immediately incorporated into any initial classical model within the 
family. A result which tells us, in particular, that quantization may be carried out more 
easily here than in a conformal quantum field theory. We end up with some final remarks. 
A brief survey can be found in f32fl. 



2 The family of spacetimes under study 

The spacetimes under study are constituted by the family of maximal spherically symmetric 
spacetimes expanded by a geodesic radial null one- form from flat spacetime (GRNSS spaces). 
They are defined through the Kerr-Schild metric 

ds 2 = dsl + 2H (r)£ ® £, (1) 

where ds 2 stands for the flat spacetime metric, H is an arbitrary function of r — the radial 
coordinate of the spherical symmetry — defined in some open region of the manifold, and 
£ is a geodesic radial null one-form. Another expression for the family is (recall that the 
signature of the metric is (—1, 1, 1, 1)) 

ds 2 = -(1 - H)dt 2 + lEdtdr + (1 + H)dr 2 + r 2 [dQ 2 + sin 2 #c^ 2 ), (2) 

where the coordinates are the customary spherical ones. In these coordinates, I = (1/ \^2) (dt+ 
dr). The other possibility, i.e. £ = (l/\^2)(dt — dr) yields the same physical results. Clearly 
enough, all the metrics have spherical symmetry and, since H = H(r), d t is an integrable 
Killing vector. In particular, for H < 1 it is timelike, for H = 1 null, and for H > 1 spacelike. 
The above coordinate choice avoids coordination problems near the possible horizons and 
Cauchy hypersurfaces, e.g. when H = 1. 

However it is worthwhile to make use of the existence of such integrable Killing vector 
in order to write down the whole family of metrics in a explicitly static form for the region 
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where H < 1. The final expression is 



ds 2 = -(1 - H)dt s 2 + — l —dr 2 + r 2 (dB 2 + sin 2 # dip 2 ), (3) 

1 — if 



where dt s is related to dt by 



cit s = dt - - -dr. (4) 

1 — ri 



This last expression for the family looks as a generalization of the Schwarzschild metric and 
may help to identify the class of spacetimes we are dealing with.Q 

The next task will be to characterize their geometrical properties. To this end we will use 
an orthonormal cobasis (local observer). In particular, we are going to choose the following, 
which also avoids any problem at the possible horizons: 

if\ , if , / if\ , if 



e 



=(l-f)d*-fdr, 0^(1 + 1)^ + 1^ (5) 

@ 2 = rd9, @ 3 = rsm6d<p. 

The Riemann tensor has the following non-zero components, together with the ones obtained 
using index symmetry (we follow the sign convention of Q]) 

7? - H " E -E - H> R -R - H> 

-ttoiOl — r - , -K0202 — -K0303 — — 7T~ , -Kl212 — -Kl313 — TT , 

if 

i?2323 = -3- (6) 

The Ricci tensor for these spacetimes has the following non-zero components 

ROO = -Rll = -\{H" + 2 ~^-\ , ^ = ^33 = 1^+^. (7) 

And the scalar curvature is given by 

^H'^JL^JL. (8) 

Finally, the Einstein tensor has for components 

G 00 = -G n = ^(V + ^), G 22 = G 33 = ~(H" + ^pJ. (9) 

Whence, we see that the Einstein tensor satisfies the relations 

Goo + G11 — 0, G22 = G33. (10) 
x t can also be written as I = (l/y/2)(dt s + dr*), where dr* = drj (1 — H) is a direct generalization of the 
usual (Wheeler) tortoise radial coordinate. 
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It is now clear that, imposing the additional relation Gn = G22, the whole set would be 
invariant under any change of the cobasis, i.e. under any Lorentz transformation (because in 
that case G a p would be proportional to g a p). In our situation, these relations are invariant 
under more restrictive conditions. Actually under any change of cobasis adapted to the 
spherical symmetry, i.e. under any Lorentz transformation between {0 ,© 1 } and between 
{@ 2 , @ 3 }. As a consequence, the physical interpretation of all these observers will be exactly 
the same. This is the reason why they are actually acceptable as consistent conditions to 
represent an admissible generalization of the absolute isotropic vacuum (G a p oc g a p)f\ — see 



e.g. [ 5UJ, 51]. Of course we could have imposed absolute isotropy in our expressions. The 
result would then have been H(r) = Ar 2 + B/r, where A, B are arbitrary constants. One 
usually sets B = in order to avoid singular behaviors near r = (i.e. the Schwarzschild 
solution), thus recovering the de Sitter spacetime for the isotropic vacuum, as expected. 
However we have preferred to keep the metric in its more general form. Actually we will 
show below how a pure de Sitter model is insufficient, for being considered as a reasonable 
model. Notice, by the way, that the spherical vacuum conditions are here obtained as a 
direct byproduct of the family of metrics under study and they do not come from any 
specific imposition. 

We shall use Einstein's field equations in the form 

= — ^r~Ta.p — Ag a /3- (12) 

Usually, one chooses units for which Gn — 1, c = 1. If the energy- matter content is of a 
quantum origin, the r.h.s. should be understood as < T a/3 >, and < A >, in accordance with 
the semiclassical approach to gravitation. 



2 For instance, in the region with H < 1, where the metric can be re-written in the explicit static form 
ds 2 = -(1 - H)dt s 2 + (1 - H)~ 1 dr 2 + r 2 (d0 2 sin 2 dtp 2 ), (11) 

we have the natural, and customarily used, static local observer associated with the cobasis 0° = 

Vl - Hdt s , 1 = (1/Vl - H)dr, 2 = rd9, 3 = r$m.6<kp, to interpret the results. 

3 The spacetimes satisfying G o + G\\ = 0, and G22 = G 33 are just of the type being described here or a 

generalization of the Nariai metric, see e.g. |}0). However the latter cannot be joined with any black hole 

model unless mass shells are introduced, which is our aim to avoid, since they introduce arbitrariness into 

the scheme. A detailed account will be given in II. 
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3 General matching conditions 



In this section we will write down the conditions for two metrics of the preceding family to 
match with each other. This kind of junction conditions constitutes by now a well-established 
branch of geometry. However, there is no canonical way of imposing them in the literature, 
and special care has to be taken in order not to confuse the reader. We will follow the 



formalism contained in p3f-|36|. 



Our focus will be in matching two metrics with the aim of having one of them as an interior 
of the other, from some point in time onwards. Here time refers to the time experienced by 
an observer inside the static region of the spacetime manifold. 

The general form of a hypersurface that clearly adjusts itself to the spherical symmetry 
of any of these spacetimes is as follows 



S : < 



'6 = 


A#, 


V = 




r = 


r(A) 


t = 


t(A), 



(13) 



where {A, X$, A^,} are the parameters of the hypersurface. Its tangent vectors are 



e\ = d\ = rd r + id t 



t: 



(14) 



where the dot means derivative with respect to A. The normal one-form is then (n ■ e*j = 0, 

i = {e,<?,\}) 



n = o(fdt — idr). 



(15) 



a is a free function in the case n is a null one- form. Otherwise, one can set n • n = ±1, by 
chosing a as 

y, ±1 



(16) 



\t 2 -r 2 + H(r -t) 2 \ 

The first junction conditions reduce to the coincidence of the first differential form of E 
at each spacetime. We must thus identify both hypersurfaces in some way. Obviously, the 
identification of (Aj)i with (Aj)2 (1 and 2 label each of the spacetimes) is the most natural 
one, due to the symmetry of the above scheme. This yields 



r(A) =0, where [/] = (/ 2 - A) s , 



r 2 



i 2 + H{r + if 



0. 



(17) 
(18) 
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The second set of junction conditions comes from (we do not use any mass shell, yet this 
could be added without problem) 



Hi 



0. 



where Hij is defined by 



„, s / d 2 <p p -pp d(pd<f>"\ 

Hin = -m o \ ^TTTTTl + J- , 



(19) 



(20) 



Here m is a vector that completes the set {e*j} to form a vectorial basis of the manifold.0 We 



will choose 



m = rdt — td r . 



This choice has the property that rh-n = a(r 2 +t 2 ). Therefore, if it vanishes, the hypersurface, 
E, becomes degenerate and the joining process itself cannot be carried out. Furthermore, 
P (A) are the parametric equations of the hypersurface ({0°, 1 , 2 , 3 } = {t, r, 9, y?}s), and 
Y p are the connection coefficients. 

Further calculations yield two supplementary conditions, namely 



[ri] 



H{tt - rr) + (1 + H)rt + (1 - H)tr - H' [f A + - + — 



f3 



r 2 t 



o. 

= 0, 



(21) 
(22) 



2 2 

where if' is the derivative of H with respect to r. 

In order to properly close the junction, it is still necessary to choose the signs of a to 
be the same in both regions, e.g. to impose a natural interior vs. exterior matching. In 
order to be sure to deal with all possible candidates, one is forced to look for all different 
solutions of hypersurfaces and spacetimes in the preceding equations. In fact, it turns out 
from Eqs. (|17|)-(f22j) that the matching conditions translate intoQ 

[r]=0, [t] = 0, [H] = 0, [H']=0. 

On the other side, H is simply a function of r, so that the two last relations may be viewed 
as a set of implicit relations defining r(A) in terms of the coefficients of Hi and H 2 . Thus, 
one arrives at the following conclusion: 

4 For the cases in which n is non-null, rh can be chosen simply as n, and Tiij becomes the second 

fundamental form. Tiij allows for dealing with a transition at the event horizon of the blak hole. 

5 There is also the possibility of having t + r = const., but this implies r = for a finite t, besides 

representing a geodesic null motion. Therefore they cannot properly descibe the final static and regular 

interior for the matter of a black hole. 
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"The only acceptable hyper surfaces fulfilling the matching conditions, that preserve the 
spherical symmetry, between two spacetimes of the GNRSS family, are those satisfying 
n(A) = r 2 (A) = R = const., h = i 2 , [H] = [H'\ = 0." 

Without losing generality, one can choose ti = t 2 = X, because of the global existence of 
the Killing vector d t . Moreover, we realize that the chosen coordinates are privileged ones, in 
which the matching is explicitly C 1 . The hypersurface, S, will be timelike, null or spacelike 
according to H < 1, H = 1 or H > 1, respectively. 

Our specific aim in the rest of this work will be to investigate if this situation can be 
actually met within our family of solutions (GNRSS space) and, also, the physical relia- 
bility of the necessary matching condition itself, including the value of H at the matching 
hypersurface. 

Of course, other solutions for the exterior picture of the object, e.g. asymptotically 
stopping, etc, may be viewed as refined treatments of the whole process. The final results, 
however, should coincide with the final (stopped) stage we will consider. 



4 The effect of the quantum vacuum. Israel's condi- 
tions. 

Ever since the formulation of the idea that a vacuum state should finish in a mode with an 
equation of state of the type p + p = 0, where p is the vacuum energy density and p its 
pressure (or stress), it seemed natural to try to match a Schwarzschild solution, acting as 
an exterior metric, with a de Sitter solution, acting as the interior one. If we had considered 
the junction conditions choosing I — n, as is often done, we would have obtained, instead 
of Eqs. (pi|) and (E3), the following ones, that account for the continuity of the second 
fundamental form, 



a{-i + H(i + f)} 



a{ ft -if + H'(f 3 - 



^fA + [3(f 2 i - i 2 f) + r 3 + t 3 ] 



0. 



0. 



(23) 
(24) 



Thus setting r\ = r 2 = R, R a constant, and H\ = 2G^m\j c 2 r ^ H 2 = (A 2 /3)r 2 (being 
mi, A 2 constants with well-known interpretations), we see that the only possibility in order 
to make them compatible is Hi = H 2 = 1. In principle this yields an interesting relation 



3 Notice that Kij = Hij taking I = n. 
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between the exterior mass and the quantum stress 



A 2 = 6GVmi/c 2 iT. (25) 

However, if Hi = H 2 = 1, it is clear that the first fundamental form becomes degenerate, 
i.e. ds| = R 2 (d\l + sin \ e 2 d\ 2 ^ ). Thus we do not have a hypersurface but actually a surface. 
This is in fact a direct consequence of the presence of a horizon (notice that we have H = 1). 
Its null character makes the above formalism insufficient, and one has to apply the junction 
conditions on a null hypersurface, or better do it once for all possibilities, as has been the 
case in the former section. This point went unnoticed in |37|], where the authors claimed 



to have matched Schwarzschild and de Sitter metrics at the horizon. In fact the true result 
eventually turns out to be negative as well, i.e. Schwarzschild and de Sitter can never be 
directly joined (this is clearly seen from the solution of the former section). This fact impelled 
other authors to consider a massive layer (usually thin) in order to overcome the problem 
§& 

A more direct demonstration of this impossibility is through Israel's conditions. Israel's 
conditions are usually associated with the continuity of the Einstein tensor in a matching 
situation. Because of Einstein's field equations, they can also be interpreted as conditions 
on the stress-energy tensor across the hypersurface. They are contained in the preceding 
ones and have proven to be useful in getting part of the set of restrictions in a very quick 
way. 

In terms of the Einstein tensor, they read 



n p G p a 



0. (26) 



In the exterior region, where we have the Schwarzschild solution, the Einstein tensor is zero, 
therefore we must have 

n P {G 2 )» a = 0, (27) 
being G 2 de Sitter's Einstein tensor which satisfies G 2 = — A 2 g, so that we get 

A 2 n = (28) 
as a necessary condition. This clearly sets A 2 = and hence a contradiction. 
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5 The case of a Schwarzschild-de Sitter exterior 



We recall that our aim was to find a regular interior of the exterior Schwarzschild geometry 
that may encompass some of the prevailing quantum ideas. However, it is not wrong to 
consider as the exterior solution a Schwarzschild-de Sitter model. That is to say, an exterior 
that may also match with recent observational evidence of a non-zero value for the cosmo- 
logical constant (see e.g. |38], as well as uncharged black holes (e.g. J5J).[]. The exterior 
metric is then generated by (Gn = c = 1) 



2mi Ai 2 

1 r , 

r 3 



(29) 



where mi is the usual Schwarzschild mass, and A x accounts for the cosmological constant 

I [3g, g. As r 



This solution is also known as the Kottler-Trefftz solution, see ET 



oo 



the spacetime tends to a de Sitter model of cosmological constant Ai (see also |20], [21] ) . 
The junction conditions translate then into 

H 2 (R) = Hl (R) = 2 -^ + ^R\ 



H' 2 {R) = H[{R) 



2mi + 9 Al R 



(30) 
(31) 



6 Regular interiors 

The preceding conditions are the ones to be imposed in order to have a proper matching of 
both spacetimes. Moreover, we will also focus on those interior solutions which are every- 
where regular. From the expressions of the Riemann tensor and the metric, we see that this 
may only be accomplished if 

H 2 (0) = 0, iZ*(0) = 0. (32) 

Thus, we finally encounter four conditions in order to have a regular interior solution. Two 
of them are imposed on the matching hypersurface, while the other two are imposed at the 
origin of the spherical symmetry. We shall thus expand H 2 at r — R or r, depending on the 
situation, or at their corresponding adimensional values r or f — 1, where f = r/R. 

7 In fact, other more general exterior spacetimes could be chosen within our family of spacetimes. For 
instance, other type of quantum vacuum contributions, see also Sect. |7.1.4| , or other classical solutions, as 
the Reissner-Nordstrom one (see Sect. |l(J). See also the Final Remarks section. 
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From now on, we will consider H 2 to be an analytic function of the variable f , a most 
natural hypothesis in view of the regular character prescribed for the interior solution.^ In 
this case, the origin conditions tell us that 

00 

H 2 (f) = £ b n r n . (33) 

71=2 

Now, one has to impose the two other conditions. Obviously it is the same to consider H 2 (r) 
or H 2 (f — 1) in the whole procedure. However, we will first work with H 2 (r — 1) in order to 
implement the junction conditions directly. From the preceding result, one immediately has 

H 2 = Y,a n (?-l) n , (34) 



n=0 



and the junction conditions tell us that 



2mi , Ai 



a = tfi(l) = -^ + y# 2 , (35) 
a 1 = ^(l) = -2(^ + ^ 2 ), (36) 

where H\{f) = (2m 1 /R)(l/f) + (A1R 2 /3)f 2 and the prime denotes now derivation with 
respect to f. 

The following step is to impose regularity of the solution — Eqs. (|32|) . We get 

00 

E(-l) B a» = (37) 

n=0 

and 

00 

J2(-l) n na n = 0, (38) 

n=0 

which, by virtue of the matching conditions, yield 

E(-l)X +2 = ^ 2 -^ (39) 

n=0 6 K 



and 



E(-l) n (^ + 2K +2 = 2 (^i? 2 - ^) . (40) 



8 It may be possible to demonstrate that analyticity or some degree of differentiability is unavoidable, 
since Hi is just a function of one variable that is only defined in a closed region, namely r G [0, R], and 
along the interval it is fixed to be finite. Otherwise, this instance serves to perform a study of the possible 
interiors. 
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It is clear that there are infinitely many possible candidates for these interiors. In the follow- 
ing sections we will analyze in more detail the properties that the big family of candidates 
share in common. 

A final remark is in order. The conditions at the origin to make the final solution "ev- 
erywhere regular", aside from the matching being also fulfilled, are important. For instance, 
H = A + Br, A and B constants, can immediately satisfy the junction conditions, therefore 
the metric is C 1 in these coordinates. However Gn and G22 yield 

Gn — — (—5 H \ G22 = , 

which is not regular at the origin, for whatever values of A or B different from zero. Of 
course, we could afford to study within these models a bigger family which also included 
singular models at the origin. In this work, as said elsewhere, we will just focus on regular 
models, because of the interest in studying their plausibility. The singular situation will be 
considered elsewhere. 



7 Isotropization 

To summarize, a basic property of any candidate to an interior solution is its behavior at 
the origin. First, due to the regularity conditions, all Riemann invariants turn out to be 
finite at the origin, as well as the components of the Einstein tensor, in accordance with 



previous works [lCj, [L^, [18], ^9[. Moreover, it is also very important to further analyze how 

they behave near the origin. 

Taking into account the expression of H2 in powers of f, we get 

1 00 

Gi^-^^ + ^ir 1 - 2 , (41) 



R 2 
1=2 



and 



G " = -wW kf '~ 2 - (42) 

It is then clear that Gn and G22 are different from each other .0 Yet we have the very relevant 
property that, for any of these spacetimes, it holds 

limGu = lim(-Goo) = limG 2 2 = limG 33 = — — . (43) 

r— »0 f— >0 r— >0 r— >0 ix 



9 It is important to notice that the conservation of the Einstein tensor implies a relation between G\\ and 
G22 of the type: G22 = Gn + G' n r/2. The same holds for the stress-energy tensor. 
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Whence, we see that a general isotropization independent of the model is actually accom- 
plished in a completely natural way. In terms of a/ we get 

o 

~M 



A 



M 



j oo 

Gn = X! Am? . 



(44) 



and 



So that 



G 



22 - — ™ 2^ o • 

M=0 Z 



(45) 



limGn = lim(— G o) = limG^ = limG 33 

f— >0 r— >0 f— >0 f— >0 

^sgC-i).^).. 



Ao 
ff 2 ' 



On the other hand, making no further assumptions on the coefficients of H 2 , we can 
isolate two of them in terms of the rest. For simplicity, we shall isolate a 2 and 03. The result 
is 



«2 



^ + ^+B-i)'y-3K 

H 3 l=i 



a 3 



6m 1 



+ £(-l)'(J-2)a,. 



1=4 



With this in hand we can write the previous expression for the central value in terms of ai, 
I > 4, 

3 



G n (0)=G 22 (0) = - 



R 2 



Jrm Ax 2 " (f- 3 )(f-2) 

IT + T^ + f J ~ 1) 2 



Z=4 



(46) 



The same can be done for the expression of if 2 in powers of r, if 2 = ^n^ r - We obtain 

3 



Gn(0) = G 22 (0) 



R 2 



it 3 /=4 



where we have used the junction conditions, i.e. 



(=2 i=2 



(47) 



(48) 



7.1 Examples 



We will consider four examples. The first two constitute the well-known proposals of JT8], [19|, 
HJ] and [^y, |21j . The other two constitute a family of new candidates that naturally arise 
from the preceding expressions. We start with the second pair 
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7.1.1 Two arbitrary powers 

We choose that only two specific powers of H(r), say M, N, be present in this case. In order 
to fulfill the regularity conditions, both must satisfy M, N > 2. However, if it holds that M, 
N > 2, the Einstein tensor, G, becomes zero at the origin. Thus, if we wish a de Sitter-like 
behavior at, and near, the origin, we must impose one of them to be equal to 2. We then 
have M = 2, N > 3 as the only suitable choice for a two power case. H 2 {f) reads 

H 2 (f) = b 2 r 2 + b N r N , N>3, (49) 



with 



»--TO^r (51) 



Gn(r) and G 22 {r) read 



Gll (f) = -A 1 + ^(^±l)(f--l), (52) 

Whence one readily sees that their finite value at the origin coincides, as expected, 

n n a Qm 1 (N+l\ 

Gn^^-Ai-— j, (54) 
and it does not change significantly if the power iV changes. 

7.1.2 Lowest powers 

This example chooses the case in which H 2 has the minimal power dependence. This amounts 
to taking iV = 3 in the preceding example or setting a\ = 0, I > 4, in the general expressions. 
Its interest lies in considering the simplest situation. The result is 

/8mi AiR 2 \ o 6mi o 
H 2 (r) = ( — + jr 2 - -^r 3 , (55) 

and, near the matching hypersurface, 
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And 



+ — ^ + J (f - l) 2 - -^(r - 1) 3 . (56) 



_ /. 24mA 24mi^ . 24m x . . . 

G n = - A x + _ + _r = -A x + -^-(r - 1), (57) 



_ / 24mA 36mi _ 12mi 36mi._ „ N 

G 22 = -(A 1 + — j+_ r =_-A 1 + — (58) 

Gn tends to — Ai as f tends to 1, in accordance with Israel's conditions. 



7.1.3 The approach of Israel and Poisson 

In reference [|18[] a plausible candidate for the energy-matter content of the interiors of regular 
black holes was proposed. Because of the impossibility of matching directly the de Sitter 
spacetime with the Schwarzschild one, the authors proposed that a singular layer of non- 
inflationary material should exist between the de Sitter core and the external Schwarzschild 
metric. The authors recognized this assumption as an attempt at overcoming the problems 
above mentioned. However the usual spirit of matching a stellar interior with a vacuum 
exterior was lost, the reason being the unavoidable presence of a singular layer, which should 
be placed close to the matching hypersurface and would act as a matter surface density. The 



authors themselves already pointed out this problem (see [44]). They considered that their 
approach could be improved by avoiding such a layer and by imposing a smooth transition 
from the hypersurface to the de Sitter core. In any case it was the only available candidate 
to continue the studies of quantum regular black holes. Further work, by these authors and 
other groups that were dealing with the same problem, addressed this idea as well as its 
physical consequences and properties (see e.g. [§, [H], 0). 

The task here will be to see whether this geometrical and physical model can be recovered 
from our analysis, solely based on the usual spirit of collapsing bodies, i.e. without involving 
any (singular) material layer. First of all, let us notice that, in our models, Goo can be 
re-expressed as 

Goo = ^(Hr)', (59) 

where ()' is the ordinary derivative with respect to r. 

We now search for a solution within our family that tends to that particular solution (see, 
for instance, [18 , |I~9| , |I~0| , |44|). This amounts to saying that we look for a de Sitter core for 
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small values of f and a quantum contribution of the type of the square of the characteristic 
curvature of Schwarzschild spacetime near the matching hypersurface, i.e. (Goo) ext. m 2 /r 6 . 
This is different to assuming that in the exterior region, close to the matching hypersurface, 
the quantum contributions turn into a cosmological-like term. We thus have a quantum 
exterior different from the one encountered in the rest of the examples. Obviously, a combined 
model could also be solved. However we prefer to maintain the spirit of the referred works 



for easier comparison of our results with theirs. In fact, at the end of Sect. |3.2| , we will show 
that the exterior vacuum contribution plays a secondary role, regardless of its precise form. 
All these features taken into account, we first put, generically, 

G °° = {B + Cr^r m 
where B and C are two constants, to be determined after imposing the matching conditions 
as well as the mentioned physical requirements. The matching conditions lead to 

' 3 H^R), 



B(B + Cr 3 ) R 2 



B{B + CR 3 ) 2 R 

where H% comes, as usual, from the external model. 

In order to reproduce the new exterior vacuum model, we see that it is possible to select 
it from our family, by setting £> ext = 0, and C~J t = ami, where a = (3Lp\, being (3 of order 
unity, and Lpi the Planck length (a is of order unity in Planckian units). f3 2 is related to 



the number and types of the quantized fields, |18|, |44j]. This choice yields 



2^ i/o^y^ 

r A r / 

where we have taken into account that the exterior region is dominated by the Schwarzschild 
geometry, for large values of f. 



Finally, B and C yield 



a mi ami 



/ 

i 



a 2 mi 

D 



(62) 



V R 3 J 



Thus, we do obtain a model within our family that copes with the desired requirements 
of references |T8| , |I9| , [Pjj . In order to check if the solution obtained coincides in fact with 



that particular model, one has to solve the problem of getting R from observed or expected 
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values of the constants involved. This will be postponed to the following section. There, the 
plausibility of all these models will be analyzed in detail. To summarize, we have proven 
here that a -everywhere smooth — spacetime model in our family satisfies all the required 
geometrical assumptions as well as the particular form of Goo corresponding to the above 
mentioned references. 



7.1.4 Dymnikova's model 

Some time after the appearance of the previous cases a new model for a regular interior of 



a black hole was proposed, see |2(| |21]| . However the approach was now different to that of 
the previous authors. In order to avoid the problem of the singular layer and to introduce a 
smooth transition between the Schwarzschild dominated region and the de Sitter central core, 
the imposition of a clear matching hypersurface was relaxed. Even more, these authors did 
not impose any matching between the two different spacetimes at all. The problem was then 
that the final picture of the collapsed object departed again from the usual, expected one. 
Now Schwarzschild was only recovered in an asymptotical sense, for f approaching infinity 
only.Q This immediately adds a natural drawback in the model, since collapsed objects have 
their matter content restricted to a finite volume. However, if a sufficiently quick convergent 
matter model can be obtained, then the lost mass, for a observer outside the horizon of the 
collapsed body could become as negligible as desired. Thus one would, at least, recover a 
trial model, interesting enough to support or reject the conclusions of the previous authors. 
In this sense Dymnikova's model is quite interesting. In a later work, together with B. 
Soltysek, they incorporated the observational fact in favor of a non-vanishing cosmological 



term in the exterior region pl| . This work includes the first pure Schwarzschild exterior 
case also straightforwardly. We shall deal in this subsection with that model, considering an 
exterior A, or Ai and a definite end of the collapsed body. 

The imposition for the energy-matter content is of the form 

Gq = G\ = 13 exp (-af 3 ) + 7, (63) 

where the constants a, (3, and 7 are found after imposing our matching conditions. First of 
10 A similar assumption is made in for the case of charged, non-rotating black holes. For the time 
being, here we will only consider those cases of major astrophysical interest, i.e. non-charged black holes. 
However, an analogous study is straightforward for their model (the core however is no longer described 
by the de Sitter spacetime or som well-known spacetime). A distinguished property of that model will be 



considered in Sect. 10 
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all we can already impose that the exterior geometry is of the Schwarzschild-de Sitter type, 
setting lim^i G*n = — Ai. This yields 



7 = — Ai — /3 exp (—a) . 



(64) 



We then integrate the expression of Gn (this is straightforward because of the assumed r 3 
dependence in this model as it was also the case in the previous model) in order to obtain 
H 2 , getting 



R 2 



a 



exp (— ar 3 ) 



+ P[/?exp(-a) + Ai 



(65) 



where we have already imposed the regularity conditions at the origin. The final step is to 
impose the matching conditions at the spatial hypersurface. They yield the condition 







exp (—a) — 1 



a 



+ exp(— a) 



Qrrii 



(66) 



Near the origin the expression of Gn tends to exp(— at)]— Ax- This constant characterizes 
the geometry of the core. Let us now assume that it is of de Sitter type as in the model of 
PH (see also next section). Therefore, it has to be equal to an internal A, say A 2 , as opposed 
to Ai, which accounts for the exterior de Sitter part. Combining both equations, one gets 

1 6mi 



exp{a) 



a 



(67) 



a[exp(a) - 1] R 3 (A 2 - Ai) 
This equation is transcendental. In order to solve it one realizes that there always exists a 
value of R, big enough, in order to make 6mi/i? 3 (A 2 — Ai) small enough. The reliability 
of this restriction of "big enough R" will be considered later. The numerical results of the 
next section will show that the associated error is indeed very small. Our aim is now to see 
whether the model considered by |20, |2l| lies also within our family of models. 

Within the restriction above, one can see that the solution is unique and appropriately 
given by 

(A 2 -A!)i? 3 



a 



Qrrii 



13 = A!- A 2 . 



(68) 



The relative error being A = a/[exp(a) — 1] ~ aexp(— a). Therefore, it can be made as 
small as needed, provided R can be taken as large as necessary. The coefficient a can be 
used to define a characteristic radius, say R c , as 



a 



6rrii 



A, 



-A, 
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3R g 
A 2 - A x 



(69) 



where R g is the Schwarzschild radius of the object. With this, for < r < R, Gu and 
G22(= — Aj + /3exp(— ctr 3 )(l — 3af 3 /2)) can be rewritten as 

G n = -Ai + (Ai - A 2 ) expi-(r/R c )% (70) 



G 22 = -Ai + (Ai - A 2 ) 



2 ViL 



exp[-(r/ii! c ) 3 ], (71) 



which exactly coincide with the expressions in the mentioned papers (see e.g. Eq. (11), (14) 
of |21|)) although they have now been obtained for the physical range < r < R. 



8 Significance and interpretation of the conditions ob- 
tained 

8.1 Physical consequences of isotropization 

Until now the precise value of the matching r (or R) has been overlooked. In part, because 
it was actually not necessary to fix it, that is, the conditions hold at any value of r, which is 
actually a rather remarkable fact. Despite the interesting observation that all the previous 
properties are independent of the precise value of R, the final aim is to match the preceding 
results with physical observations and with well-known theoretical behaviors of the quantum 
state near the origin, as well. 

Isotropization is a general property, as we have pointed out, and it can be interpreted 
as a geometrical one. By this we mean that it is actually the Einstein tensor that becomes 
isotropic, regardless of the spacetime chosen, among our general family. If we write the 
Einstein tensor in terms of the energy-momentum tensor and the A contribution, we realize 
that, at the origin, the following relation must hold 

G(0) = ^T(0)-A 2 s(0). (72) 

cr 

That is, both the vacuum contribution and the rest of the contributions of the excited states 
of the quantum field add directly. However, we can break this indeterminacy with the help 
of some physical model of quantum behavior near the origin. According to current ideas, we 
expect that the contribution from the vacuum will dominate those of the rest of states. Thus 
we simply impose T(0) to be negligible or zero when compared with Pa 2 (= (c 2 /87tGn)A2). 
On the other side, G(0) is in all cases a value that depends on R, mi, and Ai only. Therefore, 
we are faced up with a general relation of the type 

G(m 1 ,A 1 -R) = -A 2 g (73) 
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or 

G n (m 1 ,A 1 ;R) = -A 2 . (74) 

This relation will allow us to find R in terms of m\ and A 2 — Ai only. The expressions of 
Gu in more general situations have been given before. 

Finally, let us add that the model chosen here can be reinterpreted as having imposed 
an isotropization of the quantum gravitational effects in favor of the quantum vacuum. It 
is a new isotropization added to the geometrical one, which in our family has a universal 
character. 

8.2 Numerical results 

According to several recent observations |38|, |39|], in what follows we shall assume that 
Pa 1 € [10~ 10 , 2 x 1CT 8 ] erg- cm -3 . However an analysis shows that essentially the same results 
will hold even if p^ is zero, or another type of quantum vacuum contribution is assumed, as 
was the case in the Israel-Poisson model. The fundamental contribution comes directly from 
the quantum gravitational model that is imposed for the core. For the reasons given before, 
we shall consider, for the time being, an internal de Sitter-like core. Nonetheless, there 
is no present agreement about the scale at which regularization could act. This is directly 
connected with the so called " Cosmological Problem" , which is the subject of intense research 
currently.^ Our strategy will be to let this constant free and then consider the different 
results that will come out. A convenient way to handle and integrate this indeterminacy is 
to set A 2 = 10 3s Ai, being 3s the free parameter that governs the scale of renormalization. 
For instance, if s is around 40, we are then considering that regularization takes place at 
Planck scales, etc. 

8.2.1 Two arbitrary powers 

The fundamental relation is 

n The possibility of dealing with two A's has not been stressed yet and will not be discussed here. Else- 
where, we will try to perform a detailed theoretical analysis on the plausibility of having Ai — A2 of the 
order of the used numerical (Ai) and observational (A2) values. 
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mi 


s = 30 


s = 40 


s = 50 


M 


10" 9 


10 -19 


1Q -29 


1O 3 M 


10" 8 


10 -18 


1Q -28 


1O 6 M 


10' 7 


lO" 17 


lO" 27 


1O 9 M 


10" 6 


10 -16 


1Q -26 



Table 1: R in cm for various astrophysical and galactic objects and different scales of regularization 
(s = 30 corresponds to a GUT's regularization scale, s = 40 to a Planckian one, etc.). In any case 
R/LR e g is much bigger than 1 (R/L^ cg ~ 10(~ 6+s / 2 )). Therefore all of them are mainly in the 
semiclassical regime. 



whence 



R = R & </M 



I N + l 
4(iV-2) : 



(76) 



where we have put mi = Mm & , m & being the Sun mass, and R Q = ^/24m /(A 2 — Ax). This 
value only depends on the regularization scale and, in fact, is in correspondence with the 
solution for a collapsed object of one solar mass in the case of the "lowest powers" model: 
Rq G [3 x 10 21 " s ,6 x 10 2 °- s ]cm. For s = 40 we get R Q G [3 x 10- 19 ,2 x 10- 20 ]cm. Yet we 
see that the object has a quantum size very far from Planckian scales even if s is bigger. In 
general R & /Lp\ > 10 13 ! Moreover this result is valid for all since for any A^ we have that 
R G [0.6, l]R®yM. It is obvious that for any astrophysical object the final properties are 
very similar. Table |l] comprises different massive objects and regularization scales. 



8.2.2 Israel-Poisson's approach. 

We have found that the corresponding model within our family must satisfy 

B = 



6- 

In this case, B~ 2 = lim^o G 00 = A 2 , so that 

R 3 



a 

a 2 rrii 



(77) 



P 2 



6-oVA2 6 - 0jA 2 L 2 Pl 



(78) 



This model clearly depends on the coefficient f3. For instance, in order to obtain a solution, 
we must have (3 2 < 36/(AzL P {). The natural scale of regularization in this model is the 
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Planckian one since from the beginning the coefficient a was related to the Planck length. 
Obviously other regularization scales would simply change Lpi by the corresponding scale. 
Using standard values for A 2 that use a Planckian regularization scale, and that j3 2 should 
be at most of order unity [[L^, we get R ~ \f~M x 10~ 20 cm. This result is in complete 



agreement with the foregoing values, even though the models are very different, from a 
physical point of view. 

8.2.3 Dymnikova's model 

In the original model there is no frontier (hypersurface) up to which the spacetime is of 
exact Schwarzschild-de Sitter type. Thus a radius cannot be defined at which the collapsed 
object ends. However, a characteristic radius, R c , is found there, and it is proven that, for 
(r/R c ) 3 >> 1, the outside matter can be dismissed in comparison with the inner one. This 
was the main reason in order to consider that such model describes approximately an object 
of quantum size. 

In building an analogue of that situation, we decided to introduce three parameters, 
namely a, (3, and 7 in order to better mimic the model for G\\. In the rest of the models 
there are only two parameters. Now we have an extra freedom. In fact, for any R > \^2R C , 
there exists a unique solution for them, provided Ai is less than A 2 , as we expect. For the 
rest of the values, there is no solution at all. Thus we are faced with an expected freedom 
in R that can be dismissed adding a further assumption. For the time being let us see that 
there are solutions which fit in a very precise way with the previous results. 

This model was solved explicitly for the case in which a was big enough so that the 
relative error, aexp(— a), was small. Let us redefine R as R = k\^MR Ql where k is an 
adimensional constant, to be found below. Thus 

y 

Rr 



a=(^-\ =4fc 3 . (79) 



The relative error is then A = 4/c 3 exp(— 4fc 3 ). Whence, we see that choosing k of order 
unity the relative error is very small (for instance, for k = 2 it is only of 10 -11 %!). Thus 
R = k\/MR &) with k of order one, are all very good solutions of our corresponding model. 



Since there is no well defined frontier for the collapsed object in the approach of [g(J, |21 
a well-defined R does not exist. There, the assumption about the characteristic size of the 
object was ascribed to R c . In fact, R c = Rf (k\/4), so that the results in |2(], ^1|] fit well with 
ours. 
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Finally, we may add a further restriction in our model in order to eliminate the freedom in 
R. Before we proceed, it is worth remarking that the above discussion was just a check of the 
issue, if a model of the type EO, ET1] was inside our family. In fact, this is only a completion 



criteria that allows us to continue dealing with this model. It is commonly thought that no 
known quantum fields can be responsible for such a G n . 

A natural way to put an end to the indeterminacy in R is by setting, e.g., 7 = 0. With 
this assumption we have chosen an option that completely separates the quantum behavior 
in both regions. Obviously, many other choices are possible. The result of such choice is 

3 = 6m 1 ln(A 2 /A 1 ) 

A 2 -A 1 -A 1 ln(A 1 /A 2 )" 1 ; 

Now, taking into account that A 2 = 10 3s A 2 , and recalling that we expect s > 40, we obtain 

R ~ 2R c y~s. (81) 

For s = 40 this yields R ~ 7R C , which again falls within the same order as all the previous 
models. This fact leads us to believe that the main properties found hitherto are rather 
general and independent of the model that describes the collapsed object. Notice finally 
that the initial considerations that came from a direct, and wrong, matching between a de 
Sitter core and a Schwarzschild exterior yielded R = R c . So we can assert that the correct 
solution modifies the naive, and wrong, expected value by a small factor. This is quite an 
amazing conclusion. 



9 Horizons and an interpretation of the regularized 
black hole 



Looking, for instance, at Eq. (30) in [51] and comparing it with our result 



g 00 = -1 + H 2 , (82) 

we realize that substituting here our corresponding H 2 for their model, these expressions 
turn out to be very similar same, except for an overall sign due to the different signatures ( 
(+,-,-,-) instead of our (-,+,+,+)). We conclude that the same structure for the horizons 
and Cauchy hypersurfaces is obtained. In |2l| the solutions are obtained by approximations 



of the exact solution, so that these results and ours are really coincident (the relative error 
being completely negligible with respect to both exact solutions). 

24 



In general, the horizons result from the cancellation of goo- Thus we are left with a 
general set of horizons. A global study, for all the candidates encountered, has not been 
carried out yet. We could focus on examples, and try then to extract some general features 
from them, but we do not find this of primary importance.^ The main point is here, in fact, 
that the matching occurs at a radius which is substantially smaller than the Schwarzschild 
radius of the object. Therefore we will always have a typical exterior, a vacuum transition 
region extending until the matching with the object happens, and a quantum-dominated 
interior, which finally converges to a de Sitter core. In the vacuum interior region and in 
some part of the quantum object, the role of t and r are not interpreted as usual (d t changes 
its character). This is the reason for adequately treating the horizons: to see where exactly 
such changes appear. But, we can still perfectly agree in ordinary physical terms without 
requiring a general resolution of the precise radii at which horizons occur. 

Another issue is that of the topology of the solutions, if they are regular black-holes, and 
their possible extensions, i.e. "universe reborn". Its general structure can be found in E5 



where it was shown that the topology of any regular black hole should be similar to that 
of a singular-free Reissner-Nordstrom spacetime. Thus, there appears a necessary topology 
change, if one deals with their complete extension. 



10 Energy conditions 

A common point in dealing with the avoidance of singularities is to show that the required 

energy conditions in the singularity theorems (see e.g. 0) fail to be valid. 

Here we will study three local energy conditions, commonly considered in the literature: 

strong energy conditions (SEC), weak energy conditions (WEC) and null energy conditions 

(NEC). SEC are related with the formation of singularities in the collapse of an object. WEC 

are directly related with the energy density measured by an observer and NEC are useful 

in order to include some spacetimes which violate the first two, but are predicted by some 

quantum models, e.g. anti de Sitter spacetime. Although an analysis of energy conditions 

helps to understand the physics of a model, one has to be cautious ascribing to them more 

relevance than they actually have. In several systems, mainly when quantum effects play a 

fundamental role, they all may be violated with less difficulty (see e.g. the review in [ffl). 

12 With respect to the other models found here, we have got that the results are rather similar to those of 
Dymnikova and Soltysek's model plj. 
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Our results also agree with those in [29], where the authors study them in general spherically 
symmetric spacetimes. 

Let {e a }, a = 0, 1, 2, 3, be a dual vector basis of the cobasis in (^), defined by & b e a = 8 b , 
b = 0, 1, 2, 3. Any timelike vector field, V, in this manifold can be represented by 

V = A b e b , (A°) 2 = l + E(^) 2 , (83) 

i=i 

where A b are some functions. 

On the other hand, from the results of Sect. |2], the Ricci tensor is 

Ricci = R 00 (@° ® 0° - e 1 ® 1 ) + i? 22 (© 2 ® © 2 + © 3 ® ® 3 ), (84) 

where (g> is the tensor product. A similar expression is valid for the Einstein tensor. 

SEC require Rvv = RabV a V b > 0, for all V. From the expressions above, we obtain 
Rvv = Roo + (Roo + #22) P 2 ) 2 + (^ 3 ) 2 ]- Taking into account Eqs. (g) and ©, R vv = 
G22 + (Goo + C 2 2)[(^ 2 ) 2 + (-4 3 ) 2 ]- Finally, using Einstein's equations, (|TJ), and the fact that 
A 2 , A 3 are free, we get 

SEC^p + p 2 >0, p 2 >0 (85) 

where p is the energy density measured by e*o, 8np = Goo and p 2 is the tangential pressure 
(or stress) of the source, 8np 2 = G 22 - This is the usual representation of SEC. However, 
the GNRSS family allows for a different, more useful, expression. Indeed, as mentioned 
elsewhere, it is easy to show that, for a GNRSS spacetime, p 2 = — (p + rp'/2), where 
()' = d()/dr. Therefore, we can eventually write 

SEC «-> p 2 > 0, p < 0. (86) 

Following analogous steps, one finds, for WEC [Gyv > 0, for all V) 

WEC <-> p > 0, pt < 0. (87) 

In the case of NEC, V is a null vector field, V ■ V = 0, and requires the evaluation of 
R ab V a V b = G ab V a V b > 0,W. One obtains 

NEC <-> pi < 0. (88) 

Thus one sees, that a necessary condition common to all of them is that the energy profile of 
the sources be a non-increasing function. Moreover, in WEC, p > conveys the positivity 
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of energy density whereas, in SEC, P2 > refers to the positivity of tangential pressures. 
The last is in fact violated from some value of r downwards, as can be readily seen from 
P2 — —{p + rp 1 ' /2) and the fact that the source is regular at the origin. Therefore, eventually, 
the singularity is not created. 



10.1 Energy conditions for different models 



It is now easy to evaluate the fulfillment or violation of the energy conditions in the models 
presented before. It turns out that WEC, and NEC, are satisfied in all of them very easily 
for any value of r (e.g. for (anti) de Sitter core, p' = 0). One only needs to impose Ax > A 2 . 
On the other hand, as mentioned before, SEC are violated in all of them. The new result 
is that this occurs far away from the regularization scale. Indeed, a detailed analysis shows 
that SEC are violated for r < Rsec, with 



Rsec 



R 



SEC 



R. R 



SEC 



— R r 



89) 



\ 4(6i? 3 -a 2 m 1 ) / ' y 3 

where, all the quantities have been defined in Sect. 7A_ and the solutions correspond to the 
two-power model, the Israel-Poison's model and Dymnikova's model, respectively. It is clear 
from these results, that SEC are violated in the most part of the object, i.e. Rsec ~ R ( see 
also lEOin. For the evaluation of the Israel- Poison's model, we have used the same numerical 



values as in Sect. 1872. 



Let us now consider the series of models in ||22||-||24j|. In [p2J the first regular models for 
non-rotating black holes were given. We remark that all these solutions are charged black 
holes, e.g. \e\/m is of order unity, where e is the charge of the black hole. This is far away from 
astrophysical observations of black holes (see e.g. jJ0|). Nevertheless, these solutions are very 
important in this issue, since they are the first solutions to regular black holes with a clear 
interpretation in terms of a known quantum field, in this case nonlinear electrodynamics. A 
direct analysis of those works shows that, unexpectedly, all of them belong to the GNRSS 
family of spacetimes. Let us choose one of the three solutions proposed. 



For instance, the one in [p4] sets, in our notation, 

2m r / .2 



H = — 

r 



tanh 



2mr ) 



The energy density can be computed from Eqs. (|9|). We obtain 

e 2 1 1 



r 4 cosh (e 2 /2mr) e 2 s 4 y 4 cosh y~ 
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where, in the last expression we have used the same notation as the authors, that is, s = 
\e\/2m and y = 2mr/e 2 . It is easy to see that pag is always positive. However, p' AG becomes 
positive from some y downwards. In fact, a direct computation shows that the maximum 
energy density is achieved for y^ 1 tanhf^ 1 ) = 2, i.e. y ~ .4842, and that p is an increasing 
function for y < y Q . Finally, taking into account p|J that |e| must be lower than 1.05m, in 
order to have a regularized solution, we can readily show that for any r < 0.2664m, that 
is, the most part of the object, all energy conditions fail, even WEC and NEC. The rest of 
proposed models behave in an identical way.0 This adds a new (elementary) example to the 
violation of energy conditions models when quantum effects play an important role (see |46[ 
for a recent review) and shows us that energy conditions help understanding the models, but 
not necessarily should bound the search for new solutions. 



11 Quantum aspects of the models 

To find the sources for the models presented is not easy, but this will be a necessary step to 
undertake in order to show the plausibility for the avoidance of the singularities. To reach 
this goal within our scheme, one has in fact to quantize the sources. The symmetries of the 
matter-energy tensor, i.e. p + p = 0, where p is the radial pressure (or stress) — 8irp = G\\— 
and the one coming from spherical symmetry, p 2 = Pz are fundamental in order to probe 
which type of quantum fields could correspond to their sources. Taking into account the 
spherical symmetry of the spacetime and, consequently, of the Einstein tensor, one easily 
checks that all classical fields adapted to the spherical symmetry satisfy this condition (see 
for instance Sect. 3.8 of f|7|l)- This bonus was indeed expected, because spherical symmetry 
is a basic geometric symmetry in the scheme. On the other hand, the other Eq., p + p = 0, 
can be rewritten in a general covariant way, without referring to any special set of observers. 
Its expression is then 

T u = T^r = 0, (90) 

13 In [^] the authors even give a graphical representation of — g u , whence direct visualization proves our 
assertions. 
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where I is the geodesic radial null direction characteristic of the GRNSS spaces. If one 
considers a scalar field, one has 



T, 



a/3 



;i - 2O(V Q 0)(V^) + (2£ - I)^(V^)(V, 



-2^<f>V a V p <f> + I £<? 



+2 



m 2 g ai 



n<p-i[G aP + 

A2 



2(n-l) 



a/3 



(91) 



where £ is a constant representing the coupling between the scalar and the gravitational 
field, <p is the (classical) scalar field, and n is the dimension of spacetime. Whence we obtain 



2O0 2 



(92) 



a Va, and where we have used the fact that t is null and geodesic. If the field 



where ()' 

were a classical one, we should have 



;i-2O0 2 



o. 



(93) 



The first thing one notices is that there is no £ that makes Eq. (|93|) be identically satisfied. 
For instance, in the conformal case, the trace of T is the central object. Then, choosing 
£ = —(n — 2)/4(n — 1), the conformal coupling, and m — 0, one gets that the trace vanishes 
for the classical field. Therefore, Eq. ( |93|) imposes now more conditions on cf) than in the 
well-known conformal case. This is a remarkable difference with respect to conformal field 
theory. Actually Eq. (p3|) can be integrated, giving two results depending on whether £ = 1/4 
or £ 7^ 1/4. In any case, they must satisfy the field equations. 

For the GNRSS family, the free scalar field equations, (□ + m 2 + £R)4> = 0, where Q is 
the d'Alembertian operator of each spacetime, may be written as (I = 0, 1, . . .) 



-(1 + H)d 2 <p + 2Hd t d r (f> + (1 - H)d 2 .4> - d r H(d t <p - d r 4>) 



+ 



m 2 + £R 



0. 



(94) 



where R, the scalar curvature, is given by 



AH 
+ +H, 

r 



(95) 



and where a separation of variables has been done using the spherical symmetry of the 
problem. Hence, <fi in Eq. (|94| ) is a function of t and r only. We notice that R and Eq. fl94|) 
may be written more intrinsically as R = (2/6 2 )(H6 2 ), where 9 = l/2r 2 is the expansion of 
£, and 

1(1 + 1) 



2d u <p-2H<f)-2H<p+ m 2 + ^R- 



0. 



(96) 



29 



where u = (l/\/2)(r — t). It turns out eventually that the solutions of Eq. (|93|) do not 
satisfy the free scalar field equations. Furthermore, the end result is that no free classical 
field satisfies Eq. d9"3|). 

In order to see whether a quantum treatment could make Eq. fl93|) be satisfied, we shall 
begin the process of quantization of the sources for the GNRSS family (we remark that 
actually the steps are valid for any Kerr-Schild-like spacetime). First, one needs to calculate 
5. In fact, each spacetime of the GNRSS family is obtained by changing H(r), and leaving 
£ fixed. This is actually how the different interiors are obtained. Obviously, since £ is null, 
we could also choose that £ changes to a multiple of it, but this is unnecessary choosing a 
convenient definition of H. Therefore, 6g^ u is simply 

8 9llv = 2{5H)£ I1 £ U) (97) 

and, because g a/3 = r] a(3 -2H£ a £ /3 } we also have 5g a(3 = -2(5H)£ a £ f3 . Finally, g = det(g af} ) = 
det(r] a/3 + 2H£J f3 ) = det(^) = -1. 

The following step is to write down the action, S[g a p], in terms of S [r) a p] ■ By definition, 



see e.g. §6.3 in 47 



2 8S 

TaH{X) = o7— r, (98) 

where S corresponds to the classical action of matter. Performing a standard calculation, 

we get 

= S[ Va p] - J T ei (g a p)8Hd n x. (99) 

Whence, 

TuW)\ = | M • (100) 

This result may be straightforwardly generalized to include any spacetime which is in a 
Generalized Kerr-Schild (GKS) correspondence with another one, i.e. g a/ 3 = g a p + 2H£ a £p 
(see e.g. ||30|| ). Notice that the GKS form clearly takes the aspect of a problem with a 
spacetime acting as a background, not necessarily flat spacetime. 

In our case, we have Tu = (in fact, this is the case of, e.g., any Kerr-Schild metric with 
a geodesic £. Therefore, the following results are also valid, for instance, for Kerr-Newman 
metrics, or some deviations of them, which play a similar role as the GNRSS here and have 
been recently proposed, pBJ). Thus, should any classical field exist, this would mean that the 
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classical action should remain invariant under a Kerr-Schild relation. We have already shown 
that this possibility is forbidden by Eqs. (|93"D and (0) or fl96|) . Nevertheless, expression (|100|) 



tells us that T# becomes now the essential quantity in the process of quantization, in the same 
footing as in the scalar case. This result points towards the structure of a quantization 
of a GKS spacetime. Its study is, as of now, in progress and will be a matter of further 
research. It is clear that a complete solution of the problem in hand depends on this issue. 

The aim is now to test if the interiors can be the result of the quantization of some 
classical field. To that end, it is enough to focus the attention on the scalar field. We will 
therefore study the possibility that < T« >= 0. This is a natural option, due to the non- 
existence of a triad (0, £, m) satisfying Eq. fl93|) and Eq. (^) at the same time. Of course, 
this is not always a solution. Namely, given a spacetime metric, g a/ 3, its Einstein's tensor, 
G a /3, will not very often correspond to a known physical source. It is a general fact of any 
classical field equations (e.g. the monopole solution of classical electromagnetism). The same 
may happen now, even after having quantized the field. Then, there still would be a place 
for a new theory that explains it, or, eventually, it should be considered as a non-physical 
solution. It is interesting that some of the GNRSS solutions are solutions of a non-linear 
generalization of electrodynamics, see [25|, which is also connected with string/M-theory. 



Coming back to the point, we have that < T« > should vanish, while (T«)ci, the classical 
contribution, is different from zero for any spacetime belonging to the GNRSS family. 

Let us compare our present situation with that of conformal field theory and the de Sitter 
spacetime. There we have that for certain values of £ and m, say the conformal coupling 
and m = 0, the classical field verifies 

(^)ci = 0. (101) 



However, one finds after a standard calculation (n = 4, % = 1), |47], |48 



< n >= d02) 

where a 2 = — l/15a 4 , and a is the "radius" of the de Sitter spacetime. In that way, the 
de Sitter geometry, i.e. R a/ 3 = Ag a p, may be interpreted, from a source point of view, as 
the semiclassical effect of the quantum vacuum associated with the conformally coupled, 
massless, scalar field. In fact, there are several other possibilities, for other couplings and 
masses, but the effect of the anomaly is most clearly seen in the conformally coupled and 
massless case. However, T u obeys the contrary pattern, i.e. (T u ) C i ^ 0, while we pretend 
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< T a0 >-- 



(104) 



that < Tu >= 0. We shall show below that this behavior actually comes from an identical 
anomaly effect. 

For the de Sitter spacetime, one has 

T 

< T aP >= -g a/3 , (103) 

where T =< >. Therefore, since de Sitter spacetime belongs to the GNRSS family, 
{Tee)ci is different from zero, yet we have < T« >= 0. And this result has been obtained 
using the same standard formalism as that of the semiclassical approach to gravity. The 
reason is to be found in the fact that, for any local observer, we can write 

/A \ 

-A 

-A 

v -a; 

and, therefore, < >= — 4A ^ 0, while < T u >= T 00 + T u = 0. Thus, the same reasoning 
that leads from (T^)ci to < > explains, for the de Sitter spacetime, that (T«)ci is non 
zero, but < T\ > vanishes. 

Finally, we remark that all these models tend to form a de Sitter core as r/R << 1. 
Thus, the de Sitter spacetime will be a good approximation for all them, lying still in 
the semiclassical regime (see table [1]). In fact, in the models of Sect. |3.2| , the deviation is 
proportional to (r/R) 3 , and the approximation is good enough, provided that r/R << 1. On 
the other hand, the semiclassical approach makes only sense as an approximation. Therefore, 
we do not need to include higher order corrections, as of now. 

Another point of view is that T a p may be written as T a p = (Td e sitter ) a p + AT a/ 3. It 
turns out that both T^ e g itter and AT satisfy the four Wald axioms (see ^9|), for a T a p to be 
regularized. Thus, we could interpret that here we are regularizing the de Sitter's part, as 
if it were a background, (i.e. the vacuum contribution above), to which one should further 
add the contributions of the "excited" states, represented by AT a/ g. All this shows, at the 
very least, that the general scheme of quantization fits quite well within the GNRSS family. 

We recall that those above are just steps towards the confirmation of a plausible quantum 
origin for all (or some) of the solutions in the GNRSS family. To completely fix the issue, 
one should either solve Eq. (^), by finding a particular set of modes, or using Eq. (|100|) . 
Regarding the resolution of Eq. (|96|), we note that a further separation of variables between 
t and r is not possible, since the solution should be valid for the region r < R, where the 
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spacetime is non-stationary, so that usual mode decomposition is useless. Another set must 
be found, compatible with Eq. (p3|), that should be matched with some of the well-known 
exterior ones, coming from the solution of Schwarzschild or Reissner- Nordstrom spacetimes. 
Complete fixing of the exact problem seems to be a difficult task. 

Let us finally resort to a different way, less fundamental but seemingly more effective, 
in order to incorporate possible quantum corrections into the models. This may serve as a 
temporary solution, yet one that has been thoroughly used in the literature to deal with the 
effects of quantum fields in gravity. The interesting fact in our context is that, in this way, 
they are easily incorporated into any classical solution of the GNRSS family. 

Indeed, it is only necessary to compute < T 00 > to obtain < T a p >. This is because 
any GNRSS spacetime satisfies < Tq >=< T/ >, < T| >=< T| >. And, additionally, 
conservation of < T a p > yields < T| >=< Tq > +(r/2) < Tq > (see also footnote |J). Thus, 
the stress-energy regularization is similar as for trivially scalar fields, where the trace of < 
T a p > is the only piece to be regularized necessarily. This greatly simplifies the computation. 
Moreover, another very interesting (and important) property of the GRNSS spaces is that the 
solution of T = T + T 2 is given by H = Hi + H 2 , where Hi, and H 2 are particular solutions 
for Ti, and T 2 , respectively. This last property can be interpreted as a manifestation of 
the equivalence principle, since this is the same as asserting that m = mi + 7712, where 
rrii are the mass function in each case. Therefore, in order to find the (well-established) 
quantum corrections, it is only necessary to integrate G[H q ] = quantum corrections and 
write H reg = H dass + H q in order to know < T a p >. 

In fact, if one just considers (as is commonly assumed, see e.g. ^7, IE, 0), that semi- 
classical corrections are, roughly speaking, of the form oc TZ 2 , with TZ 2 some geometrical 
Riemannian invariant of the classical spacetime, then the corrections are obtained through 
direct integration of 

^(H q r) >r = XU 2 , 

where A is of order GTi, i.e. of order unity in Planckian units, and is related with the number 
and types of quantized fields, [IB, <0[. An application of this method to two well known cases, 
namely de Sitter and Schwarzschild, will speak for itself. In the de Sitter case we have that 
R = 4A = const, therefore R 2 — as well as any other Riemannian invariant — is a constant 
and the quantum correction of a macroscopic de Sitter metric turns out to be a de Sitter 
term. Thus one arrives at H rcn = A ren r 2 /3 (where the subscript "ren" means renormalized) . 
This is a well known result for the renormalization of a de Sitter spacetime (see e.g. H7], pjfl). 
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The Schwarzschild case has R = 0, for it is a classical vacuum solution. However, R a p\^R a ^ ^ 
is 48m 2 /r 6 . Thus, we set a quantum vacuum contribution of the type Am 2 /r 6 . Integration 
yields the remarkable result that the quantum vacuum corrections due to its polarization 
yield a spacetime metric belonging to the GRNSS spaces with H = 2m/ r — Am 2 /r 4 . Other 
sources are to be found within the effective actions of M-theory or string theory (see e.g. 
p2| and references thereinF]) and could be easily incorporated into our scheme, once the 
corresponding Einsteinian metric is computed. We will consider them elsewhere. 



12 Final remarks 

The first thing to be noticed is the intrinsic freedom present in our model, which is as large as 
the measure of the set of analytic functions of one variable. This is a very rewarding feature, 
since it allows to impose further restrictions coming from new proposals. In particular, 
it surely will be a helpful tool when trying to find explicitly a quantum field responsible 
for the Gn and G22 in the funadamental uncharged case. For comparison, in all previous 
works, based on a unique model, the prospective of finding a quantum field related with 
their energy-matter content was hopeless. 

The second thing to be mentioned is, that we have here been able to develop in depth 
some relevant issues left incomplete in other works. For instance, in [[0|, the authors studied 



the interpretation of the de Sitter region as a new universe birth. They also considered the 
effects of the evaporation of the black hole which would be interesting to be computed for 
the models presented here. In ||44||, the authors studied the stability of the model. In the 



present paper, all these statements have been substantially implemented, and this in a very 
natural way (see the preceding sections). 

On the other hand, the source origin for both examples in the section before is clear, but 
they do not provide the right solution of the problem we have in hand. The de Sitter one 
is not enough, since a de Sitter interior cannot be directly joined with an exterior coming 
from a black hole. And the Schwarzschild solution is singular at the origin. Then, one 
has to find a deviation of the de Sitter spacetime (a non-isotropic metric) still physically 
meaningful, although this metric will be only needed macroscopically. This will be the 
subject of subsequent work (Refs. p], O, contain hints towards this direction). 



We have here given just the first steps towards the goal of the (semi-classical) quantization 
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Indeed, it is not difficult to see that all those solutions belong to the GNRSS class. 
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of any suitable source. Our preliminary results show undoubtedly that this process, for the 
case of Kerr-Schild type of metrics, is under the same footing as when one has to deal with 
a conformal field theory. 

We have also shown that the solutions behave as constituting an anomaly effect of their 
sources. Moreover, we have seen how quantum field corrections can be easily incorporated 
into any initial classical model within the family. A result which tells us that the quantization 
scheme may be carried out more easily here than in conformal field theory. All these findings 
not only compel us to believe that black hole singularities are likely to be removed by 
quantum effects but, in our view, open also a new window for the search of a compatible 
quantum field that, once regularized, may yield the same result for, at least, a particular T 
inside these models. 

Finally, the preceding sections were centered on the free field equations. The possibility 
of considering some potential is certainly important, if it can be related to some physical 
theory, as in |28| for the case of rotating black holes. We reckon that the field approach 
above may yield a different and useful strategy towards the same goal. 
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